Phase dynamics and particle production in preheating 
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Abstract 

We study a simple model of a massive inflaton field <j> coupled to another scalar filed x with 



interaction term g 2 (f> 2 x 2 - We use the theory developed by Kofman et at. for the first stage of 
preheating to give a full description of the dynamics of the x field modes, including the behaviour 
of the phase, in terms of the iteration of a simple family of circle maps. The parameters of this 
family of maps are a function of time when expansion of the universe is taken into account. With 
^vqj ' this more detailed description, we obtain a systematic study of the efficiency of particle production 

as a function of the inflaton field and coupling parameters, and we find that for g < 3 x 10~ 4 the 
£f) • broad resonance ceases during the first stage of preheating. 

> 

m 

w> ■ 1 Introduction 

O , 

The success of the reheating stage after inflation is crucial to most realizations of the inflationary 
paradigm. The universe needs to recover a temperature high enough for primordial nuclesynthesis to 
take place in accordance to the usual pattern of the standard cosmological model. Excluding scenarios 
designed to avoid the extreme cooling produced by inflation (such as, for instance, the warm inflation 
scenario), it is important that the necessary post-inflationary reheating be efficiently achieved. 

The reheating mechanism was proposed as a period, immediately after inflation, during which the 
inflaton field 4> oscillates coherently about its ground state and swiftly transfers its energy into ultra- 
relativistic matter and radiation, here modelled by another scalar field X- This processes depends on the 
coupling between (f> and x m the interaction Lagrangian. The classical theory of reheating was developed 
in El ID H] • The importance of broad resonance and approximations to deal with non-perturbative 
effects were introduced in IH1 El GUI ■ The theory was put forward in El E21 i where the analysis 
included the effects of expansion of the universe, see also JHlEj- This represented a shift away from 
the simple picture of static Mathieu resonant bands, due to large phase fluctuations, which behave 
irregularly in the non-perturbative regime in an expanding universe. Extensions in this setting include, 
among others, the study of metric perturbations and the application to string or supersymmetric theories 

[T3II1II3QJQJI2S]. 

The present understanding of the process (see [21]) distinguishes two parts in it. A preheating 
mechanism by which fluctuations of the inflaton couple to one (or more [22]) scalar fields, inducing the 
resonant amplification of perturbations in the latter. Depending on the coupling, efficient energy transfer 
requires the amplitude of the inflaton oscillation to be rather large, away from the narrow resonance 
regime where only perturbations with wave numbers in small intervals are unstable. As the amplitude 
of the perturbed field grows, back reaction effects may have to be considered, since the frequency of 
the inflaton oscillations is no longer given by its mass, but depends also on the total number density 
of the perturbed field particles through the coupling term. The first stage of preheating is the period 
when these back reaction effects are negligible, and the inflaton field dynamics is approximated by its 
uncoupled equations. Preheating ends when resonant amplification terminates, either because of the 
decreasing amplitude of the inflaton field or due to the back reaction and rescattering effects of the 
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second stage. After the second stage of preheating, the reheating period corresponds to the decay of 
the perturbed fields as well as that of the inflaton field that comes out of the preheating period, leaving 
the universe after thermalization with the temperature required by the subsequent processes, namely 
nuclesynthesis. 

We consider a basic model describing the inflaton field <j> interacting with a scalar field x in a flat 
FRW universe 

L = + ^ X ,iX A - V{4>) - Vintifr x). (1) 

This is the simplest model that still contains the basic features for the understanding of particle creation 
in the early universe and one of the few models for which an analytical study can be performed, see 
also (HlElEBj- With this in mind we concentrate on the simplest chaotic model with the potential 
V{(f)) = l/2m^(j) 2 and interaction potential Vi nt ((f>,x) — 9 2 4> 2 X 2 - The evolution of the flat FRW universe 
is given by 

3H ' 2 = 5; + v ® + ¥ 2 + 5 w ) (2) 

where H = R/R and R is the FRW scalar factor. The equations of motion in a FRW universe for a 
homogeneous scalar field <fi coupled to the fc-mode of the x field are given by 

4> + 3Hi>+(ml + g 2 xt)^ = (3) 
Xk+ZHxk+uftt)xk = 0, (4) 

where u>l(t) — k 2 /R 2 + g 2 cj) 2 . 

The rate of production of particles of a given momentum k is determined by the evolution of the 
perturbed field mode Xk- The number density nj~{t) of particles with momentum k can be evaluated as 
the energy of that mode divided by the energy of each particle 

nm = ^(\Xk\ 2 ^ 2 (t)\Xk\ 2 )-\. (5) 

The exponential growth Xk(t) oc e Alfc * leads to exponential growth of the occupation numbers nk{t) and 
the total number density of x-particles is given by 

n x (t) = j^ Jdkk 2 n k (t). (6) 

The problem of determining the efficiency of particle production for a given model is thus reduced to the 
evaluation of /i^ as a function of the parameters of the model. In general this is not an easy task, and 
most estimates are based on numerical integrations for typical parameter values 

However, the theory developed in [T| yields analytic results in Minkowski space time, which provide 
an approximate simplified model that works rather well in the FRW scenario. The starting point of 
the theory is the fact that, for the model <j2j, 10} , preheating requires the amplitude of the inflaton 
field's oscillations, to verify g$(t) > m<f,. This means preheating is dominated by broad resonance, 
and the theory is based on the approximations that hold when g$(t)/m0 3> 1. The conclusions hold 
independently of the detailed form of the inflaton potential away from its minimum. 

In this paper we extend the formalism of to give a full description of the dynamics of the phase 
of the field modes Xk, which in turn determines the evolution of the growth factor \x k . In Section [2 we 
consider the case of a non expanding universe and show that the behaviour of the phase can be described 
in terms of the iteration of a simple family of circle maps. The orbits of this family are of two possible 
types, depending on the value of the perturbation amplitude, and their asymptotic behaviour, which 
determines the growth factor, is always independent of the initial condition. 

In Section UH we show that the results that hold for Minkowski space time can be used to model the 
phase and growth factor evolution in a FRW universe. Expansion may be considered simply by taking 
the parameters of the Minkowski equations to be prescribed functions of time. Using this model as an 
alternative to numerical integration of the full equations, we check the estimates given in the literature 
for the total number density of x particles created during preheating and for the duration of the first 
stage of preheating as a function of the interaction parameter g. 
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2 Phase dynamics in Minkowski space-time 



The important differences between the narrow and broad parametric resonance mechanisms in the context 
of cosmological models with post-inflationary reheating have first been noticed in the analysis of the 
evolution of the amplitudes of the Xk modes in a static universe 0E1ISI- Although the equations in 
Minkowski space time lack some of the fundamental ingredients to understand the overall efficiency of 
the reheating process, they can and indeed they should be considered as a toy model that sheds light on 
the mechanisms at play when the expansion of the universe and other effects, such as back reaction and 
rescattering, are taken into account. 

This was the approach followed in Kofman et al. 1^, where an analytic theory of broad resonance 
in preheating was established, relying on a detailed study of broad parametric resonance driven by 
the harmonic oscillations of an inflaton field without expansion of the universe. This study is directed 
towards the computation of the k-mode growth factors and the phase dynamics of the Xk modes is 
not explicitly derived. However, Kofman et al. do mention the basic features of this phase dynamics, 
and they use them to explain the characteristics of mode amplification in an expanding universe, which 
they dubbed 'stochastic resonance' in order to stress the difference with respect to the usual resonant 
bands scenario of the Mathieu equation. 

In this section we explicitly compute the phase evolution equations in Minkowski space time and show 
that phase stochasticity is already present in this model, as one of the two possible dynamical regimes, 
together with the fixed phase behaviour identified in pp. In this random phase regime, in the complement 
of the resonant bands, the growth factor [ik is effectively zero, but a typical orbit undergoes random 
sequences of amplitude amplifications and reductions, much like in the case of stochastic resonance. This 
phenomenon is pointed out frequently in the literature [251 1261 l2"T] and can be described analytically as 
one of the consequences of the theory when we consider the global phase dynamics of the Xk modes QQ[£]. 

We shall start by recalling the method of Kofman et al. to approximate in the broad resonance regime 
the solution of equation <Q} in Minkowski space-time 



Xk + ul(t)xk = 0, 



(7) 



where Lj\{t) = ak + b sin 2 (t) and the time variable is now t — > m^t. The parameters au and b are given by 
a-k = k 2 /m 2 ^ and b = g 2 A 2 /m^, where A is the constant amplitude of the field 4>. Typical values of the 
parameters are g 2 < 10~ 6 , m = 10 -6 m p 
In broad resonance, Vb^> I. 
Let the Xk{t) be of the form 



A = 



where < a < 1, and thus b < a 2 x 10 6 d EU E 



Xk(t) 



Oik 



exp — i 



w fc (s)ds 



0k 



V^Mt) 



where ak and (3k are constants. Introducing |JH| in to (JJ we obtain 



exp i / w fe (s)ds 



(8) 



Xk +u 2 k {t) 



1 -i d 1 
4 I k dt 



1 d 

2uik dt 



i d 1 
-77 m LU k 

dt 



Xk = 0. 



So (SJ approximates the solution of Q provided that 

-id 



-id / , -id. 

W,, — \ LO h — hi LUk 

k dt \ k dt 



« 1, 
< 1 



(9) 

(10) 
(11) 



hold. These are the adiabatic conditions identified in |2!)| . In the present case, \Jb ^> 1 and the adiabatic 
conditions are fulfilled except in the neighbourhood of tj — jir, j = 0, 1, . . ., when 4>(tj) — that is, 
every time the inflaton field crosses zero. Hence, the breakdown of the approximation given by JHJ occurs 
periodically, and an approximate global solution of Q can be constructed from a sequence of adiabatic 
solutions 



xl(t;alPl) 



exp 



Lo k (s)ds 



Pi 



exp i / cj fc (s)ds 



(12) 
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where the parameters {pel, p k ) for consecutive j will be determined by the behaviour of he solution of 
for t close to tj. In a small neighbourhood of tj = jir, j = 0, 1, . . . equation can be approximated by 



Xk + (a k + b{t - t 3 ) 2 ) Xk = 0. 



(13) 



Equation ll 1-311 has exact solution in the form of parabolic cylinder functions |30| . The asymptotic 
behaviour for large t of this solutions is of the form JSJ, and this provides the relation between the 
coefficients (ak,Pk) of this adiabatic approximation on either side of tj. Following Kofman et al. pp, 
this relation can be written as 



■j+i 



- 1 


R* K i 








1 


-~D~k 


~D* 



(14) 



where 9 k — f Q 3 w(s)ds. The complex numbers R K and D K are given by 1/D K = ^/l + p\ exp(i^5 K ) and 
R K /D K = -ip K with p K = exp(-7TK 2 /2), k 2 = a k j\fb, and 

'1 



¥>« = arg T 



-' 1 + ln ^ 



(15) 



The parameter n = k/ \J Agm^ £ [0, 1] is the normalised wave vector of the mode, and ■J Agm^ the cut- 
off wave vector introduced in pQ . Since the number density of \k particles with momentum k is equal to 



|/3fe| 2 , one can use ltT4l and ltl5ll to calculate the number density of particles 



tj in terms of n k 



/4 



J+ 1 



\(3 k \ 2 . The growth index p J K , defined by n k +l = n J k exp(27T/z£), is given by [Q 
— In (l + 2p 2 - 2p K y/l + f% M~<Pn + arg^ - argaj[ + 29{) 



or, in terms of the phase u k 



arg/# 



9 3 k of the field Xk when t = tj, 



2tt 



In (l + 2p\ - 2p K ^l + p%sm(-<p K + 2^)) 



after 



(16) 



(17) 



In FigureQ]we show, for b — 10 3 and a k — 1, the analytic curve for the growth factor p\ as a function 
of the phase v given by II17H and the numerical curve p K — p K (v) obtained from the numerical integration 
of the full equation along half a period of the inflaton field. We see that, as v varies in [0, n], p takes 
positive and negative values. The phase interval for which p is negative depends slightly on the value of 

K. 

As hinted by Kofman et al. 1,, equations (11411 and 11511 can be used to obtain the dynamics for the 
phase iP k . In the remaining part of this section we explicitly compute the map v J k +1 — Pb, K {v k ) and show 
that it can be approximated by a simple family of circle maps. From 11411 one gets 

ii +1 = e{b,K) 



arg 



ple-^e^ 



(18) 



where 9{b, k) 
given by 



,J+ 1 - 



f Q uj(s)ds. An approximate expression for the phase map up to terms of order n 2 is 



s/2 sin \A — cos v 



V/2 log^- 




Ci cos V, 



C2 cos v k sin v k 



c 3 sin vi 



4v / 2cos isL svuvl 



(19) 



where a = 2.074 - log 2/k 2 , c 2 = -1.363 + 1.4141og2/K 2 , and c 3 = -0.147 - log2/« 2 . 

The maps 1181 . II19II for b = 10 3 and k = 0.5 are shown in Figure [2 together with the phase map 
obtained from the numerical integration of the full equations for the same values of the parameters 
over half a period of the inflaton field. 
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The properties of the family II18H and its approximation l|19fl are best understood by looking at the 
behaviour of the family Pb.o{v) parametrised by v&, 



. . rr v2sini^ — cosi^ 

P M (v)=2Vb + atan : (20) 

V2 cos v — sin v 

This family of circle maps is periodic with period 7r in the parameter Vb, and its bifurcation diagram 
for Vb e [0, 7r] is shown in Figure^! We see that the map has two different regimes. For tan 2\/b G [—1,1], 
the map has a strongly attractive fixed point, and the phase converges rapidly to this fixed point (typically 
in a couple of iterates). For the remaining parameter values, the phase orbit has random oscillations 
around the mean value that varies between tt/8 and 7r/2 — 7r/8, or between n/8 + ir and 3ir/2 — n/8. The 
fixed point equation 

cos 2v „ r- , . 

= tan2V6 (21) 



V2 — sin 2v 



is satisfied for two values of v for each Vb S [— 7r/8,7r/8]U[7r/2 — 7r/8, 7r/2+7r/8], one of which corresponds 
to an unstable fixed point and the other to the stable fixed point shown in the bifurcation diagram of 
Figure The derivative of P^o at the stable fixed point is 

KM = - l. ~ < i. (22) 

3 — 2v2sm2z/ 

where the stable fixed point v 6 [— 7r/2 — 7r/8,7r/8] U [tt/2 — 7r/8,7r + tt/8]. The equality Pl (v) = 1 is 
obtained at the boundaries of the random phase region, but the derivative decreases rather sharply into 
the stable fixed point region, where P b ' Q (T>) < 0.5 for most values. 

The asymptotic value of growth factor [iq for the k = mode as a function of b can then be computed 
from equation ifTTjl evaluated at the fixed point or averaged over the random orbits, for values of Vb e 
[— 7r/8, 7r/8] U [tt/2 — 7r/8, tt/2 + tt/8] or in the complement of this interval, respectively. In Figure 31 
we show a plot of the asymptotic value of as a function of b computed analytically from equations 
(1201) . i!17|l . as described, and numerically from the integration of the full equations 0. We see that 
the stability regions identified in P correspond to the random wandering of the phase v of the field \ 
at consecutive tj — jir around its average values. The random phase regimes correspond to different 
probability measures in the phase interval, for all of which the average value of (i is zero. 

For other values of k, the global dynamics shares the qualitative properties of the family Pb,o- In 
Figure we show the same information of Figures and H obtained for P M with 300 < \/b < 304 and 

K = l/2. 

Equation lfl8ll provides a good approximation to the exact phase dynamics even for moderate values 
of b. In Figure El we show the approximate phase map II18J1 and the phase map computed numerically 
from the full equation Q for b — 10 and K = 1/2. This, together with the rapid relaxation rates of 
the phase dynamics in the fixed point regime, is the reason why equations Ill7|l . I|18ll are still useful to 
obtain estimates for the growth number in FRW space-time. The existence of two simple regimes for the 
phase dynamics, one of them characterised by rapid relaxation to a fixed point and the other by random 
wandering of the phase with a well defined probability density, shows that the system keeps no record of 
the initial phase and, in this sense, has no memory. This explains why the evolution of the occupation 
number is independent of the initial phase, while the growth factor per period y? K depends strongly on 
the phase v 3 k . We shall come back to this point later. 



3 Dynamics in an expanding universe 

In this section we will show that equations l|17H , (I18H can be used to compute the growth factor of the x 
field modes during the first stage of the preheating period in a flat FRW. 

It is well known pQ, that the coherence of the oscillations of the field is not disturbed until the 
energy density of the x field significantly contributes to Q , © and back- reaction and rescattering effects 
start to change the mechanism of the growth of the occupation number of the produced particles. More 
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precisely, Kofman et al. show in pQ that these effects are negligible during the first preheating stage, 
that ends when the total number density n x of \ particles satisfies 

n x (t) « (23) 

where is the varying amplitude of the inflaton field 4>. 

In the first stage of preheating, equations EU decouple from <@J|, and the evolution of the inflaton 
field and of the scale factor R(t) is given in good approximation by 

m = maat , *(t) = J« t ) = (| ) 2/3 . (24) 

Equation can be reduced to the form through the change of variable X k = R 3 ^ 2 Xk, yielding 

X k +w&t)X k = 0, (25) 



where vo 1 — — 9 — + '■ ~ 1 5-<5 and the last term is very small after inflation and will be 



g 2 Htf , i 

disregarded. 

We see that the \k modes are still governed by an equation of the form like the one we considered in 
Minkowski space time, but now the parameters a k and b change with time. By definition, the preheating 
period ends when g$(t)/m^ ~ 1, and so, during preheating, the rate of variation of those parameters and 
the oscillations of the inflaton field are much slower than the oscillations of the Xfc modes. As pointed 
out in [I] , the basic assumptions for the approximation developed for Minkowski space time are thus still 
valid in preheating, and the changes in occupation numbers n k will occur at t = jir with exponential 
growth rate given by 11711 , provided that the decreasing amplitude of the perturbations and the redshift 
of the wave numbers are taken into account. Hence, Kofman et al. model particle production in the first 
stage of preheating through 



< = ^ In (l + 2p% - 2p Kj + p2 . sin(-^ 3 + 2^)) , (26) 



where Kj = k/(R(tj) y j / gm ( f > ^(tj)). 

We may also think of the phase dynamics as being essentially governed by the iteration of equations 
11811 . but now the parameter b decreases in time, crossing the strips associated to the random phase and 
fixed point regimes, and giving rise to non trivial phase dynamics. Kofman et al. describe the effects of 
the variation of the inflaton amplitude as implying random phase dynamics, and build their estimates 
for the total number density evolution n x (t) from Il21ill . treating the phase v 3 K as a random variable. 

We shall extend the approach of [J to study the phase dynamics in an expanding universe, taking 
the phase iteration map to be given by 



4 +1 



B{bi,Ki) + arg (^p^e-i^ - :/<.. e - ) , (27) 



where ^fb] — g^it^lm^. 

Equations ffijl , (l2fill , provide an alternative to numerical integrations of the full equations do compute 
the occupation number of a given mode as a function of time. The phase, growth factor and occupation 
number for a typical orbit (we have taken bo — 5 x 10 3 and kq — 0.1) as obtained from the iteration of 
equations 12711 . J 2 fill are shown in Figure where the values given by numerical integration of equations 
(1251 for the same initial conditions and parameter values are also plotted. Iteration and integration 



were carried out until preheating ends with \/b(tj) sa 1. We see 'reminiscences' of the phase dynamics 
of the Minkowski model. In particular, the fixed point regime interval is clearly visible after the first 
few <f> oscillations. Also shown are the values of these same quantities averaged over the initial phase 
v k . We see that, due to the characteristics of the phase dynamics, the efficiency of particle production 
is insensitive to the initial phase of the field mode Xfei m spite of the strong dependence of the growth 
factor per period on the phase v 3 k . 
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We shall now use equations 112 711 , 12 (ill above to look at the total number density n x (t) and check the 
estimates given in PQ. We consider the contributions to n x (t) of every mode such that k € [0, fc*(0) = 
TJgW(ff)m$] at t = 0, and compute for each j the leading mode's growth factor, /J, , given by 



[ w°e[0,27r] J 



with 6^ ; = 5 x 10 3 ($(<j)/$(0)) 2 , k 2 = k 2 $(0)/ ($(i,)-R 2 (t,)), and the j-th iterate of the map p7|. 

where the parameters k and 6 must be updated at each iteration. Then, whereas the estimate in pQ is 



(ffm $(O)) 3 / 2 
647r 2 i? 3 (i)V0.137ri 



n x (t) = _ yy Zt7. irrri ^ ex P( 2 x °- 13 ^ ( 29 ) 



we have 



„(«,). ( -""'^;° ))3/2 exp 2„£ ff! , (30, 



or 



%fe)= (g^jog exp ( 2 ^„A (31) 



647r 2 i? 3 (t)^/7r 2 ^ =1 ^ V i=i 

if instead of actually determining the leading mode we assume that it corresponds to /c*/2. The curves 
corresponding to equations Q2S)fl. ll-'UH . ll-'illi are shown in Figure |H1 We can still improve on the estimates 
given by ll3f)|l or l(3"T)l , Since 

- 14^4) I (^)d/v, (32) 

we may use (29 j EEt to compute this integral numerically from rik(tj) — 1/2 exp fJ- l Kj ^j ■ The 

result is also shown in Figure|Hl and we see that for this value of g the estimates of pQ are very accurate. 

Finally, we have used again 12711 , 12611 together with i-VIl to obtain a systematic study both of the 
behaviour of n x (t) as a function of the parameter g, and of the duration of the first stage of preheating 
as defined by ll23ll . as a function of g. The results for g in the range 10~ 4 < g < 10 _1 are shown in 
Figure El where we have also plotted the result of the estimates of for these quantities. We see that 
the duration of the first stage of preheating depends rather irregularly on the parameter g, and that 
the estimates of pQ (full line in the figure) provide a good lower bound for most values of g. The least 
squares linear fit (slashed line) yields a slightly larger value for the typical duration of the first stage. 
On the other hand, the value of g below which preheating always ceases during the first stage was found 
to be g = 3 x 10~ 4 , as predicted in [[]. 



4 Conclusion 

In the simplest preheating scenario, where the coherent oscillations of the uncoupled inflaton field drive 
the amplification of the mode amplitudes of a field %, we consider first the broad resonance regime in 
Minkowski space-time and use the theory of scattering in parabolic potentials developed in pQ to obtain 
the map whose iteration governs the phase dynamics of the modes Xfe • It is well known that the phase 
dynamics, the consecutive values of the phase of the Xfc fields at the times tj when the inflaton amplitude 
goes through zero, determine the growth rates of the modes. In this work we show that the features of 
this phase dynamics are given by the properties of a simple family of circle maps. The orbits of this family 
of maps are of two types, rapid convergence to fixed point solutions, and random oscillations around 
an average value. Hence, the 'stochastic resonance' identified in [[] in the dynamics of an expanding 
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universe is also present in the absence of expansion. The fixed phase and stochastic regimes occur in 
consecutive intervals of the value of the forcing amplitude. In the first case, the fixed point is always 
associated with a positive value of the growth factor fij = l/{2iz) log (nP +1 /vP) that controls the growth 
of the number of particles -n? in each mode for t = tj. Thus, in this case, the average growth of the 
occupation numbers of the modes is exponential. In the second case, we show that the phase sampling 
is always such that the average growth factor is zero. 

We then consider the case of an expanding universe, with the assumptions that hold in the first stage 
of preheating, and show that the equations for the phase dynamics and the growth number derived for 
Minkowski space time still provide a good approximation of the true solutions, once the decay of the 
inflaton amplitude is taken into account. Moreover, the qualitative behaviour of the phase and growth 
number evolution is reminiscent of the behaviour found in the case without expansion, in the sense that it 
can be interpreted as a random phase regime followed by a slowly varying phase regime where occupation 
number growth is approximately exponential. These two regimes occur as the inflaton decay slows down 
and the perturbation amplitude crosses more and more slowly the intervals that give rise to fixed phase 
behaviour. 

We use this approximation to obtain a systematic study of the behaviour of the total number density 
of created particles over time, and of the end of the first stage of preheating as a function of the 4> — x 
coupling parameter g. Comparison with the estimates presented in show an overall good agreement. 
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Figure 1: For b = 10 3 and au = 1 (hence k 2 = 10 3 / 2 ) the analytic curve for the growth factor as a 

function of the phase given hv pmiatinn HI 7| l fsnlirl lineV and the numerical curve nhtained from the 

integration of the equal 
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Figure 2: For b — 10 3 and k — 0.5 the maps l|18H (solid line) and II19|I (dashed line). Also shown is the 
phase map obtained from the numerical integration of the full equations Q (full circle) for the same 
values of the parameters. 




Figure 3: Bifurcation diagram of the family of circle maps ll2?Hl for \fb g [0, n]. 
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Figure 4: The asymptotic value of fio as a function of b for \/b £ [10tt, lln] computed analytically from 
equations p0|. itTTIl (full line) and numerically from the integration of the full equations (0 (dotted 
line). We have taken j ma x — 200, and the two lines almost overlap. Also shown (in grey) are all the 
values of/4, j = 100, 101,... ,200. 
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Figure 5: (a) Bifurcation diagram for the map P M with Vb e [300, 304] and k = 1/2. (b) The asymptotic 
value of fj, K as a function of Vb for the same values of the parameters computed analytically from equations 
i!18|l . 11711 (full line), and numerically from the integration of the full equations Q (dotted line). We 
have taken j max — 200, and the two lines almost overlap. Also shown (in grey) are all the values of [x J K , 
j = 100, 101,..., 200. 
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Figure 6: The approximate phase map (solid line) and the phase map computed numerically from equa- 
tions Q (dotted line) for b = 10 and k = 1/2. 
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Figure 7: For b = 5 x 10 3 and ko = 0.1, the phase (a), growth factor (b) and occupation number (c) 
with initial conditions corresponding to n® = 1/2 and — 0. The values obtained from the iteration of 
equations 1270 . 12fifl are plotted as full circles, and the values given by numerical integration of equations 
(l25l) for the same initial conditions and parameter values are plotted as open circles. Iteration and 
integration were carried out until the end of preheating when y/b(tj) « 1. Also shown are the values of 
these same quantities averaged over the initial phase (full triangles for the iterated maps QJJ, lf2fijl 
and open triangles for the numerical values). 
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Figure 8: For bg = 5 x 10 3 (hence g — 0.3 x 10 -3 ) the curves corresponding to equations H29I1 (full line), 
(13011 (dotted line), l-Tlll (slash-and-dot line) and Il-'i2t (slashed line). The evolution was computed until 
the end of preheating when v^>30 — 1- 




-log 10 g 

Figure 9: Duration of the first stage of preheating j\ as a function of g for 10~ 4 < g < 10 _1 (dots). The 
linear least square approximation (slashed line), are also shown the estimate given in (full line) and 
the cut-off curve defined by bj = 1 (dotted line) . The values of j\ that lie above the cut-off curve bj = 1 
were not considered in the least square approximation 
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